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Detecting  Intersection  Among  Star  Polygons 

Delfin  Y.  Montuno 
Alain  Fbumier 
University  of  Toronto 


ABSTRACT 

We  give  algorithms  to  solve  in  0(n)  preprocessing  time  and 
0(n)  space  the  following  problems  involving  star  polygons,  where  n 
is  the  number  of  vertices:  Given  a  point  and  a  star  polygon,  find  in 
log(n)  time  where  the  point  is  with  respect  to  the  star  polygon. 
Given  a  line  and  a  star  polygon,  find  in  O(logn)  time  an  intersect¬ 
ing  edge  of  the  star  polygon  if  they  intersect  or  the  neeirest  vertex 
or  vertices,  otherwise.  Given  a  line  segment  and  a  star  polygon, 
find  in  O(logn)  lime  a  point  of  common  intersection  if  they  inter¬ 
sect  or  a  locally  nearest  vertex  or  vertices,  otherwise.  Given  a 
convex  polygon  and  a  star  polygon,  find  in  O(logn)  time  a  point  of 
common  intersection  or  a  locally  nearest  vertex-pair,  otherwise. 
Given  two  star  polygons,  find  in  0(n)  lime  a  poinl  of  common 
intersection  if  they  intersect  or  the  nearest  neighbor  vertex-pair, 
otherwise. 
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1.  Introduction 

Numerous  results  have  been  published  in  computational  geometry  in  solv¬ 
ing  convex  intersection  problem  [Sha76,  Sha78,  Mul78].  A  variation  of  this 
problem  is  the  detection  problem  [Cha80,  DobBl].  Sublinear  time  algo¬ 
rithms  have  been  found  for  problems  involving  convex  objects  in  2-  and  3- 
dimensions.  Some  attempts  at  solving  similar  problems  involving  non- 
convex  objects  can  be  found  in  [Sha78].  The  usual  strategy  for  solving 
problems  involving  a  convex  object  and  a  non-convex  object  is  to  decom¬ 
pose  the  non-convex  object  into  its  convex  parts  and  then  to  perform  on 
them  the  intersection  or  detection  computation.  Very  often,  the  problem 
becomes  a  problem  of  decomposition  [Avi80,  Cha79,  Fen75.  Kei82, 
Sac8l.Sch78].  The  solution  to  such  a  decomposition,  however,  is  not  neces¬ 
sarily  an  optimal  solution  to  the  intersection  or  detection  problem.  For 
example,  in  the  detection  of  intersection  between  a  star  polygon  S  and  a 
convex  polygon  C  (see  definitions  below),  both  with  0(n)  vertices,  we  would 
have  to  decompose  S  into  its  convex  parts  and  in  the  worst  case  we  could 
obtain  0(n)  of  them.  Thus  a  detection  algorithm  with  O(logn)  time  com¬ 
plexity  for  two  convex  polygons  with  a  total  of  0(n)  vertices  would  require 
O(nlogn)  time  to  solve  our  detection  problem.  One  of  the  reasons  is 
because  the  structural  property  of  the  non-convex  object,  in  the  present 
case  a  star  polygon,  is  not  fully  exploited. 

The  purpose  of  this  paper  is  an  attempt  to  fill  this  gap  in  presenting  results 
on  2-dimensional  detection  problems  involving  a  non-convex  object  and  a 
convex  object,  where  the  non-convex  object  belongs  to  the  class  of  star 
polygons. 

Using  0(n)  storage  and  0(n)  preprocessing  time  for  star  polygons  with  0(n) 
vertices,  we  are  able  to  detect:  1)  in  O(logn)  if  a  point  is  inside  a  star 
polygon  or  not;  2)  in  O(logn)  time  if  a  line  intersects  a  star  polygon  or  not; 
3)  in  0{logn)  time  if  a  line  segment  intersects  a  star  polygon  or  not;  4)  in 
O(logn)  time  if  a  convex  polygon  intersects  a  star  polygon  or  not;  and  5)  in 
0(n)  time  if  two  steir  polygons  are  intersecting  or  not.  Additional  informa¬ 
tion  such  as  nearest  vertex  to  an  edge  and  nearest  neighbor  vertex-pair  of 
two  star  polygons  are  also  found  for  non-intersecting  situations. 

The  first  problem  has  been  solved  by  [Sha78],  but  we  present  it  here  with 
minor  modifications  for  the  sake  of  completeness.  The  corresponding 
intersection  detection  problems  involving  convex  polygons  have  also  been 
solved  [Cha80,  DobSl].  We  note  here  that  computing  the  intersection  of  two 
star  polygons  takes  O(n^)  time  in  the  worst  case,  where  n  is  the  total 
number  of  vertices  [Sha78]. 
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We  use  the  definitions  and  observaiions  in  Section  2  and  describe  the  5 
detection  algorithms  in  Section  3.  In  Section  4.  we  describe  how  0(n) 
preprocessing  time  and  0(n)  storage  space  can  be  obtained  for  fast  detec¬ 
tion  algorithms.  In  Section  5.  we  analyze  and  prove  the  complexity  of  the 
detection  algorithms.  Finally,  we  conclude  with  some  further  problems 
and  research  directions. 

2.  Definitions  and  Observations 

2.1.  Definitions 

Def.  1 

The  convex  hull  vertices  of  a  simple  polygon  is  the  subset  of  the  vertices  of 
the  polygon  that  appear  as  vertices  of  its  convex  hull, 

Def.  2 

Two  points  in  a  simple  polygon  are  mutually  visible  iff  the  straight  line  con¬ 
necting  the  two  points  is  completely  inside  the  polygon. 

Def.  3 

A  kernel  of  a  simple  polygon  is  a  convex  polygon  which  is  the  possibly 
empty  subset  of  points  of  the  polygon  such  that  all  points  in  the  polygon 
are  visible  from  all  points  in  the  kernel. 

Def.  4 

A  kernel  point  of  a  simple  polygon  is  a  point  in  the  kernel  of  the  polygon. 

Def.  5 

A  star  polygon  is  a  simple  polygon  with  a  nonempty  kernel.  (Henceforth 
abbreviated  as  star  and  represented  as  a  sequence  of  counter-clockwise 
ordered  vertices,  (^'o,  Vi,  ....  where  n  is  the  number  of  vertices.)  See 

Fig  1  for  an  example  of  a  star  polygon,  its  kernel,  kernel  point,  and  convex 
hull  vertices. 

Def.  6 

Two  geometrical  objects  intersect  iff  one  or  more  points  belong  to  both  of 
them. 

Def.  7 

A  swe2p(k,  Vi,  Vj)  of  a  star  with  respect  to  kernel  point  k  is  the  region 
swept  by  rotating  about  k  line[k,t/i)  to  line[k,  Vj)  including  the  lines  as 
boundary,  where  Iine[k,7;)  is  the  half-infinite  line  starting  at  k  and  contain¬ 
ing  V :  and  i  <  j  in  counter-clockwise  order. 

Def.  B 

A  span(k,  Vj)  of  a  star  is  a  sweep(k.'i;^,i»j)  of  the  star  such  that  the  two 
vertices,  and  Vy,  are  consecutive  convex  hull  vertices  of  the  star. 
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Pig.  1  Star  polygon  S.  its  kernel  K(S),  kernel  point  k,  and  convex  hull  ver¬ 
tices 


Def.  9 

A  sector(k,  v^,  v^)  of  a  star  is  a  sweep(k,i>^,uy)  of  the  star  such  that  the  two 
vertices,  and  Vj,  are  consecutive  vertices  of  the  star.  See  Fig.  2  for  an 
example  of  sweep,  span,  and  sector  of  a  star  polygon. 

Def.  10 

A  polygonal  chain  pc(i;i . vj)  of  a  simple  polygon  is  a  connected  subse¬ 

quence  of  vertices  of  the  polygon  from  Vi  to  vj. 

Def.  11 

A  span  chain  of  a  star  is  a  polygonal  chain  of  the  star  contained  in  a  span. 
For  example,  a  span  chain  of  the  star  polygon  in  Pig.  2  is  ...,  ug), 

Def.  12 

A  convex  chain  of  a  polygon  (polygonal  chain)  is  a  polygonal  chain  of  the 
convex  hull  of  the  polygon  (polygonal  chain).  For  example,  a  convex  chain 
of  span(k,U3,^l9)  in  Tig.  2  is  (i^g,  Vg). 

Def.  13 

A  clockwise  convex  chain  of  a  span(k,i;i  ,Uj  )  is  any  one  of  the  convex  chains 
ending  at  from  Vi,  where  i  <  i’  ^  j.  (Note  that  though  the  convexity  of 
the  chains  is  either  clockwise  or  counter-clockwise,  the  chains  are  always 
listed  in  counter-clockwise  order.)  Counter-clockwise  convex  chains  of  a 
span(k,x'i,i/j)  are  defined  similarly.  For  example  in  span(k.vg,ug)  of  the  star 
polygon  in  Pig.  2,  (^'s.'i's.'u?)  and  (vq.v^.Vq)  are  some  of  the  clockwise  convex 
chadns  and  {vQ,V'i,vg)  and  (v^.Ds.vg)  are  some  of  the  counter-clockwise  con¬ 
vex  chains. 
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DeL  14 

The  closest  clockwise  convex  cliain  of  a  star  to  an  outside  point  p  is  the 
clockwise  convex  chain  of  the  star  that  intersects  line[p,k]  closest  to  p, 
where  k  is  a  kernel  point  of  the  star.  The  closest  counter-clockwise  convex 
chain  of  a  star  to  an  outside  point  p  is  defined  similarly.  The  closest 
counter-clockwise  convex  chains  of  the  star  in  Fig.  2  to  points  pi  and  p2 
are  {vQ.vj.vg)  and  {vs.Vg),  respectively. 


f  ;  VO 
CHt,  Va, 


Fig.  2  Sector(k.Vi,'y2).  span(k,T;3,VQ),  sweep(k.i'2.'t'9).  and  convex  chains  of  a 
star  polygon. 


Def.  15 

A  visible  chain  (ui,  V2 . Vp)  of  a  convex  polygon  C  with  respect  to  a  point 

k  is  a  polygonal  chain  of  C  such  that  all  points  on  the  chain  are  visible 
from  k.  If  k  is  in  C.  then  it  is  C.  See  Fig.  3. 

Def.  16 

The  nearest  vertex  of  a  visible  convex  chain  to  a  reference  point  and  a 
reference  line  is  called  the  vertex  of  partition.  If  the  chain  and  the  line 
intersects,  then  the  vertex  of  partition  is  the  farthest  point  on  the  chain 
from  the  line  on  the  side  containing  the  reference  point.  Otherwise,  it  is 
the  nearest  point  on  the  chain  to  the  line.  The  head  (tail)  portion  of  the 
visible  chain  is  the  part  of  the  chain  that  is  before  (after)  the  vertex  of  par¬ 
tition  in  the  counter-clockwise  order.  See  Fig  3.  Note  that  there  could  be 
at  most  two  such  vertices  i/i+ij.  In  such  case  the  convex  chain  is  parti¬ 
tioned  into  3  parts:  head  portion,  (  ....  Vi)-.  tail  portion,  (i^i+i,  ...);  and  line 
segment  portion  (vt,  'L'i+i). 


-5- 


/ 


Tig  3 
tion  V 


4- 


2.2.  Observations 
Obs.  1 

The  convex  hull  of  a  star  polygon  with  n  vertices  can  be  found  in  0(n)  time 
[BhaSl]. 

Obs.  2 

The  kernel  of  a  simple  polygon  can  be  found  in  0(n)  time  [Lee79]. 

Obs.  3 

Star  polygon  can  be  trtcingulated  with  one  steiner  point  in  the  kernel  in 
0(n)  time. 

Obs.  4 

Point,  line  segment,  and  convex  polygon  belong  to  a  special  class  of  star 
polygons  for  which  their  kernels  are  the  polygons  themselves. 

Obs.  5 

A  line  intersects  a  star  polygon  iff  it  intersects  its  convex  hull. 

Obs.  6 

A  line  segment  intersects  a  star  polygon  iff  one  of  its  endpoints  is  inside 
the  star  polygon  or  it  intersects  an  edge  of  the  star  polygon. 
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Obs.  7 

A  convex  polygon  intersects  a  star  polygon  iff  one  of  its  vertices  is  in  the 
star  polygon,  one  of  the  star  polygon's  vertices  is  in  the  convex  polygon  or 
one  of  its  edges  intersects  some  edge  of  the  star  polygon. 

Obs.  6 

Two  star  polygons  with  n  vertices  each  could  have  at  most  0(n*)  intersec¬ 
tions  [Sha78],  while  a  star  polygon  and  a  convex  polygon  with  n  vertices 
each  could  have  at  most  0(n)  intersections.  Similarly,  a  star  polygon  with 
n  vertices  and  a  line  could  have  at  most  0(n)  intersections. 
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3.  Algc»iliiDi8 

We  now  present  the  5  intersection  detection  algorithms.  In  practice,  the 
preprocessing  step  is  done  once  and  for  all  but  it  is  included  here  with  the 
detection  step  so  as  to  indicate  the  kind  of  preprocessing  needed  for  each 
algorithm.  So  for  all  the  algorithms,  both  the  preprocessing  and  detection 
steps  will  be  given  together.  They  will  be  evaluated  separately  in  Sections 
4  and  5. 

3.1.  A1  Detecting  inclusion  of  a  point  in  a  star  polygon. 

INPUT: 

A  point  p  and  a  star  polygon  S,  Vi .  where  vertices  are 

arranged  in  counter-clockwise  order  and  n  is  the  number  of  vertices. 

OUTPUT; 

YES.  if  p  is  in  S  and  NO,  otherwise.  In  either  case,  indicate  the  sector 
where  p  is. 


PREPROCESSING: 

1.  Rnd  kernel  of  S,  K(S). 

2.  Choose  a  point  k  in  K(S). 

3.  Store  the  vertices  of  S  in  a  structure 
suitable  for  binary  search  with 
respect  to  k. 

DETECTION: 

1.  Determine  by  binary  search  which  sector  p  is. 
Point  p  is  in  sector(k,Ui.i»i+i)  iff  the  oriented 
emgles  (p,k,i»i)  and  (p.k.i^t+i)  have  opposite  signs. 
(In  this  way,  we  do  not  have  to  compute  any 
polar  angles.) 

2.  Output  sector(k,i;i,i;i+i)  where  p  is. 

Ifp  is  inside  triangle(k,i;i,T;i4.i)  then 

output  YES  else  NO, 
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3.2.  A2  Detectiqg  mteraecUtiB  jbetinLUii  a  tine  cffid  a  star  polygon. 

INPUT: 

A  line  I  and  a  star  polygon  S,  (sq.  . ^n-i).  where  the  vertices  are 

arranged  in  counter-clockwise  order  and  n  is  the  number  of  vertices, 

OUTPUT: 

YES  and  an  intersecting  edge  of  S  if  intersecting;  NO  and  the  nearest  ver¬ 
tex  or  vertices,  otherwise. 

PREPROCESSING: 

1.  Find  convex  hull  of  S,  CH(S),  (hq,  . 

where  m  <  n. 

2.  Find  kernel  of  S,  K(S). 

3.  Choose  an  interior  point  k  of  K(S). 

4.  Store  the  vertices  of  S  in  a  structure  suitable 
for  binary  search  with  respect  to  k. 

5.  Store  the  vertices  of  CH(S)  in  a  structure  suitable 
for  binary  search  with  respect  to  k. 
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mBdlON: 

1.  Determine  by  bin«a:y  search,  the  span(k,i/y,Vj+i) 
containing  the  perpendicular  from  k  to  Z . 

2.  If  convex  hull  edge{vj  ,vj  ^.i)  intersects  I  then 

Output  YES  and  an  edge  of  S. 

i 

else 

\  /•  Find  limits  of  binary  search  according  to 
where  the  line(i'j  .Uj  +  i)  intersects  I  with 
respect  to  the  perpendicular  from  k  to  f .  ♦/ 
/*  below  and  above  refers  to  before  and 
after  in  counter-clockwise  order.  */ 
if  line(vy.x'y+i)  intersects  I  below  (above) 
the  perpendicular  or  I  is  in  CH(S) 
and  parallel  to  edge{vj  then 


Let  sector(k.^<i,^'i+.l)  be  the  sector  containing 
the  line  parallel  to  f  at  k  and  below  (above) 
the  perpendicular  from  k  to  Z . 

Let  the  search  limits  be  {vi,Vj)  for  below, 
for  above,  and  either  of  the  two 
if  I  is  parallel  to  edge(i'j-,i;j  ^-i). 

Binary  search  between  vertices  {vi.Vj) 

(or  for  above)  until  either  an 

edge  (i;j' .i/i-4.i)  intersects  I  orvi-  =  +  i 
and  Vi'  is  not  on  I . 

If  edge  (vf  .Vf  4-1)  intersects  I  then 

CXitput  YES  and  an  edge  of  S  or  a  vertex  of  S  if 
Vi’  is  equalled  to  + 

else 

Output  NO  and  the  nearest  vertex  or 
vertices  to  I . 

i 

I 

else  /♦  edge(i;^-,i;y+i)  is  parallel  to  L  and 
is  outside  of  CH(S)  */ 

f 

Output  NO  and  the  nearest  vertices  to  I . 

! 

i 
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3^.  A3  Detectmg  interseeticm  between  aline  segfeent  and  a  star  polj^gen. 

INPUT: 

A  line  segment  Is  and  a  star  polygon  (i;o.  where  the  vertices 

are  arranged  in  counter-clock\^se  order  and  n  is  the  number  of  vertices. 

OUTPUT: 

YES  and  the  sector  where  one  of  the  endpoints  Is  in  or  an  intersecting  edge 
of  S,  if  they  intersect.  NO  and  the  nearest  vertex  (vertices)  in  the  sweep 
where  Is  is  visible  from  a  kernel  point  of  S,  otherwise. 

PREPROCESSING: 

1.  Find  convex  hull  of  S,  CH((S). 

2.  Find  kernel  of  S.  K(S). 

3.  Choose  a  point  k  in  K(S). 

4.  Order  the  vertices  of  S  with  respect  to  k  in  a 
structure  that  allows  binary  search. 

5.  For  each  span(k,■L'i.^'J-).  find  the  set  of 
clockwise  convex  chains  and  the  set  of 
counter-clockwise  convex  chains. 
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DEXlSnCOi: 

1.  if  either  endpoints  of  Is  is  in  S  then 

I 

Output  YES  and  the  sector  where 
the  endpoint  of  ts  is. 

i 

2.  else 

2.1  Let  sector and 
sector(k,^'i,^»i^.l)  be  the  sectors  containing 
the  endpoints,  lep  and  rep,  of  Is ,  where 

j  >  i  and  (k,lep,rep)  is  a  clockwise 
oriented  triangle. 

/*  convex  chcdns  are  used  to  facilitate  binary  search 
on  a  sequence  of  relevant  edges  of  the  star  polygon 
associated  with  vertices  on  the  chains.  */ 

/*  intersection  means  intersection  only  with 
the  original  edges  of  S.  */ 

2.2  if  Is  intersects  clockwise  convex  chain  from  vx  or 

Is  intersects  counter-clockwise  convex  chain  from  vp, 
where  vx  is  Vj  if  line[rep,lep)  intersects 
line[i;j  andi^^+i  otherwise 
andi»p  is  if  line[lep,  rep)  intersects 
line[t;,+i,T;i)  and  Vi  otherwise  then 

Output  YES  and  one  of  the  intersecting  edges. 

i 

else 

2.3  if  sectors  containing  lep  and  rep  are  in 

different  spans  then 
use  a  modified  version  of  algorithm  A2. 
else 

I 

Output  NO  and  a  vertex  (vertices)  on  one 
of  the  two  convex  chains  nearest  to  Is . 

i 

i 

i 


See  Hg.  4. 
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Rg  4  Intersection  detection  of  line  segments  and  a  star  polygon. 


3.4.  A4  Detecting  intersection  between  a  convex  polygon  and  a  star  polygon. 
INPUT: 

A  convex  polygon  C,  (cq,  Cj .  c^_i)  and  a  star  polygon  S  (sq,  Sj . 

where  both  sets  of  vertices  are  arranged  in  counter-clockwise  order  eind  m 
and  n  are  the  number  of  vertices  in  C  and  S,  respectively. 

OUTPUT: 

YES  and  a  pair  of  intersecting  edges  or  a  common  point,  if  they  intersect. 
NO  and  the  nearest  vertex  of  S  to  C  in  the  'sweep’  where  C  is  visible,  other¬ 
wise. 


Same  as  for  Algorithm  A3. 
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i^racnoN:  / 

1.  if  So  (or  k)  is  in  C  then 

CXitputTES  and  Sq  (or  k). 

2.  else  if  cq  is  In  S  then 

I 

Output  YES  and  cq. 

i 

3.  else  if  C  intersects  CH(S)  and  contains  a  vertex  of  S  then 

Output  YES  and  the  enclosed  vertex  of  S. 

4.  else 

4.1  Find  the  visible  chain  VC,  (1,  ....  r)  of  C 
with  respect  to  k  in  span(k,X,/9), 

WLOG,  let  VC  =  (1,  Ci,  Ci4.i.  ....  r). 

/*  The  case  in  which  there  are  2  vertices  of 
partition  can  be  solved  in  the  same  manner 
shown  below  with  additional  intersection 
detection  for  the  line  segment  portion  of 
the  visible  convex  chain.  */ 

4.2  WLOG,  assume  there  is  only  one  vertex  of  partition. 
Partition  VC  into 

Ic  =  (1,  Ci,  . c*)  and 

rc  =  (c|j.  . r)  with 

reference  point  k  and  line[X,p], 
where  is  the  vertex  of  partition. 

4.3  Let  sector {k,sy,Sj+j)  be  the  sector 

containing  c;^. 

4.4  if  is  in  triangle(k,Sj  ,Sj+i)  then 

Output  YES  and  • 

i 

4.5  else 

I 

/■*■  Convex  chains  are  used  to  facilitate  binary 
search  on  a  sequence  of  relevant  edges  of  the 
star  polygon  associated  "with  vertices  of  the 
convex  chains.  */ 

/*  intersection  means  intersection  only  with 
the  original  edges  of  S  */ 

4.5.1  if  rc  intersects  clockwise  convex  chain  from  v  or 

Ic  intersects  counter-clockwise  convex  chain  from  v , 

where  v  is  if  line[c^;^-l.c;b) 

intersects  lin€[sj,sy4.i)  ands^+i  otherwise  then 
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Output  YES  and  a  pair  of  intersecting  edges. 

i 

4.5.2  else 

I 

Output  NO  and  the  vertex  or  vertices 
nearest  VC  from  one  of  the  two  convex 
chedns. 

i 

i 

i 

See  fig.  5. 


Fig.  5  Intersection  detection  of  a  visible  chain  from  point  k  and  a  star  po¬ 
lygon. 
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3.5.  A5  Detect!^;  intersectkm  betweraa  tiro  star  poljgons. 

WPUT: 

Two  star  polygons  P,  (po.Pi . Pn-i)  and  Q,  (go.gi .  where  both 

sets  of  vertices  are  arranged  in  counter-cloclcwise  order  and  n  and  m  are 
their  respective  numbers  of  vertices. 

OUTPUT: 

YES  and  a  pair  of  intersecting  edges  or  a  common  point,  if  they  intersect. 
NO  and  the  nearest  neighbor  vertex-pair,  otherwise. 

PREPROCESSING: 

For  both  P  and  Q  do  the  preprocessing  steps  1  to  4  of  A3. 

DETECTION: 

1.  if  po  is  in  Q  then 

Output  YES  and  Po- 

I 

2.  else  if  go  is  in  P  then 

1 

Output  YES  and  go- 

i 

3.  else 

I 

Find  the  'visible  chain’  VP  (VQ)  with 
respect  to  A:^,  a  point  in  K(Q) 

(kp ,  a  point  in  K(P)). 

while  (scanning  VP  in  counter-clockwise  order 
and  VQ  in  clockwise  order  until  end  of 
either  chain  is  reached  or  an  intersection 
is  found) 

I 

determine  the  nearest  neighbor  vertex  pair. 

I 

if  intersection  is  found  then 

t 

Output  YES  and  a  pair  of  intersecting  edges  or 
a  common  point. 

I 

else 

I 

Output  NO  and  the  nearest  neighbor  vertex-pair 
in  the  scan  order. 

i 

i 
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4.  l¥e|HT>cesaiog  time  and  Mxxrage  spaee 

The  preprocessing  steps  of  a  star  include:  finding  its  kernel,  finding  its  con¬ 
vex  hull,  storing  its  vertices  in  a  data  structure  suitaWe  for  binary  search, 
and  finding  counter-clockwise  and  clockwise  convex  chains  in  all  its  spans. 

From  observations  1  and  2,  the  preprocessing  steps  of  algorithms  A1  and  A2 
require  only  0(n)  time  and  0(n)  space,  where  n  is  the  number  of  vertices  of 
a  star.  However,  it  remains  to  show  how  the  convex  chains  can  be  found  in 
0(n)  time  and  stored  in  0(n)  space.  Note  that  in  the  worst  case,  there 
could  be  0(n)  convex  chains  each  of  length  0(n).  Thus  we  can  not  explicitly 
determine  and  store  them  and  yet  expect  to  obtain  0(n)  space  and  time 
complexity  (see  Fig.  6). 


list  cWkiivi^^  CTfi/Hit 


tV,,  Va,  U4,  ) 

t  Ui  ,  Ui ,  V4  ,  Ufc  ) 
iv. ,  .  V4,  Uff ) 

C\j,,  Vi,  04) 

C  Vi ,  '^4,1/5) 

CV,,  \}^) 


fig.  6  An  example  of  a  star  polygon  containing  convex  chains  of  O(ti^) 
length  if  they  are  all  explicitly  listed. 


The  strategy  we  will  follow  is  to  explicitly  list  each  vertex  exactly  once  and 
store  only  branching  information  about  overlapping  convex  chains.  (See 
Rg.  7.)  We  only  give  here  an  algorithm  for  determining  the  clockwise  con¬ 
vex  chains  in  a  span  of  a  star  polygon  since  each  span  chain  can  be  pro¬ 
cessed  similarly  and  independently  of  the  others  (and  the  counter- 
clockvnse  convex  chains  in  a  given  span  are  found  similarly). 
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4.1.  Ag  <stedkwi»e  4eoTOrex  rhMrein  a  apan  of  a  star  poiyj|<ML 

wpxrr. 

A  span  chain  SC  =  (i;i .  Vj^)  Ln  counter-clockwise  order  and  a  kernel 

point  k  of  a  star  polygon  S. 

OUTPUT: 

The  set  of  clockwise  convex  chains  in  SC  stored  in  ao'rays  INFO[l:m]. 
LOWER_LlMlT[l;m].  and  BRANCH_FLAG[l:m]. 

NOTAnONS: 

INFO[i]  is  a  pointer  to  where  the  information  about  vertex  is  stored  in 
linear  arrays  LOWER_LIMIT[*]  &  BRANCH_FLAG[*]. 

BRANCH_FLAG[INFO[i]]  indicates  how  vertex  Vi,  among  other  vertices 
which  are  also  visible  to  Vj-,  is  related  to  the  vertex  Vj,  where  Vj  precedes 
Vi  in  counter-clockwise  order  in  the  span  containing  both  of  them.  It  is  a 
non-branch  point  (0)  if  angle(k,i;j  ,i/i)  is  maximum  among  all  vertices  visi¬ 
ble  to  Vj  and  a  branch  point  (l),  otherwise.  This  means  that  if  a  vertex  Vi 
is  a  branch  point  then  it  is  the  start  of  a  non-first  convex  chain  which 
branches  out  from  vj.  This  also  implies  that  there  is  at  least  one  convex 
chain  with  a  sequence  of  non-branch  points  from  Vj  and  maybe  some  other 
convex  chains  each  starting  with  a  different  branch  point. 

LOWER_LlMIT[INFO[i]]  indicates  the  lower  limit  of  the  convex  chain  ending 
at  Vi.  In  other  words,  it  is  the  location  of  the  nearest  branch  point  (vertex) 
to  Vi  from  to  1  with  respect  to  the  convex  chain  containing  it. 
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MEfHOD: 

1 .  If  m  <  1  then 

if  m  =  0  then  Output  no  chain, 
else 

I 

INFO[l]  =  1: 

L0WER_UMIT[1]  =  1; 

BRANCH_FLAG[1]  =  1; 
size[l]  =  1; 

i 

Stop. 

i 

2.  /*  empty  stack  and  initialize  ♦/ 
empty  stack  SK; 

/*  there  is  one  temporary  list  for  each  vertex 
containing  pointers  to  the  vertices  visible 
from  it  and  it  is  initially  empty  */ 
append  Vg  to  the  list  of  vertices  for  Uj, 

/*  size[i]  stores  the  number  of  vj's  visible  from  ♦/ 
size[l]  =  1; 

push  V 1  and  then  ug  into  SK, 

/*  j  is  the  index  of  the  vertex  being  processed  ♦/ 
j  =  3: 

3.  size[2:m]  =  0; 
while  (j  ^  m) 

\/*  append  all  vj's  to  the  list  of  vertices 
for  Vi  if  angle (k.Vj.Vj)  is  maximum  for  all 

Vj,  ...,  Vi_i,  Vi+i . Vj_i  for  1  <  i  <  j  <  m  ♦/ 

Let  Vi  be  the  top  element  of  SK. 

Let  Vg  be  the  second  element  from  the  top  of  SK. 
while  (SK  is  not  empty  and  has 
more  than  one  element  & 
angle (k.Vj, Vi)  <  angle  (k.Vj.Vj,)) 

pop  SK; 

Update  indices  i  and  s ; 

j 

append  Vj  to  the  list  of  vertices  for  vi . 
size[i]  =  size[i]  +  1; 
push  Vj  into  SK; 

j  =  j  +  1: 

j 

4.  /*  Store  convex  chains  in  INF0[l:m],  LOWER_LIMIT[l:m]. 

andBRANCH-FLAG[l:m].  */ 
empty  stack  SK; 
pos  =  1; 

BRANCH  =  1; 

NON_BRANCH  =  0; 

INFO[l]  =pos; 

LOWERJJMIT[pos]  =  1; 
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BRANCH_FLAG[pos]  =  BRANCH; 
push  Vi  into  SK; 
ajRRENT_BRANCH  =  pos; 

5.  Let  Vi  be  the  top  element  of  SK. 

’w^ile  (SK  is  not  empty  8c 

the  list  of  vertices  for  is  not  empty) 

Let  Vj  be  the  last  vertex  in 
the  list  of  vertices  for  Vi. 
pos  =  pos  +  1; 

INFO[j]  =  pos; 

if  the  number  of  vertices  pointed  by  Vi  =  Bize[i]  then 

i 

/♦  first  vertex  removed  from  from  list  of  */ 

LOWERJLIMIT[pos]  =  CURRENT_BRANCH: 

BRANCR_FLAG[pos]  =  NON_BRANCH; 

i 

else 

I 

LOWER_LIMIT[pos]  =  INFO[i]: 

BRANCHJPLAG[pos]  =  BRANCH; 

CURRENT-BRANCH  =  pos; 

i 

push  Vj  into  SK; 

delete  vj  from  the  list  of  vertices  of  v^. 
while  (SK  is  not  empty  &:  the  list  of  the 
top  vertex  is  empty) 

pop  SK; 

i 

Let  Vi  be  the  top  vertex  in  SK. 

i 

Steps  1,  2,  and  4  take  a  constant  number  of  operations,  while  steps  3 
and  5  take  0(m)  time  since  each  vertex  is  processed  a  constant 
number  of  time.  (Note  that  for  each  j  in  the  loop  only  vertices  in  the 
stack  are  processed,  and  that  at  most  one  vertex  is  pushed  down  the 
stack.)  This  implies  that  the  processing  for  counter-clockwise  and 
clockwise  convex  chedns  in  all  spans  of  a  given  star  polygon  with  n  ver¬ 
tices  can  be  done  in  0(n)  time  and  store  in  0(n)  space.  See  Rg.  7  and 
Table  1  for  an  example  of  0(n)  storage  of  clockwise  convex  chains. 

TTius.  all  preprocessings  can  be  done  in  0(n)  time  and  require  0(n) 
space. 
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¥ig.  7  0(n)  storage  structure  of  convex  chains. 


5.  Detection  cocoplexity 

It  is  straightforward  to  show  that  O(logn)  time  is  sufficient  to  execute 
the  detection  step  of  algorithms  A1  and  A2,  where  n  is  the  number  of 
vertices  of  the  star  polygon.  As  for  algorithms  A3  and  A4.  their  time 
complexity  are  respectively  O(logn)  and  0(log(n+m)),  where  n  amd  m 
are  the  numbers  of  vertices  of  the  star  polygon  and  the  convex 
polygon,  respectively. 

In  the  worst  case,  algorithm  A3  requires  6logn  +  O(logn)  +  k  time  — 
2logn  for  step  1;  21ogn  for  step  2.1  (to  determine  the  sectors  contain¬ 
ing  the  endpoints  of  a  line  segment):  8logn  for  step  2.2  (to  find  the 
intersection  of  Is  with  convex  chains):  and  O(logn)  for  step  2.3  — 
where  k  is  some  constant.  Thus  the  detection  step  of  algorithm  A3 
requires  O(logn)  time. 

In  the  worst  case,  Algorithm  A4  requires  logn  each  for  steps  2  and  4.3: 
logm  for  step  4.2  and  2logm  each  for  steps  1  and  4.1:  0(log(n+m))  for 
steps  3  and  4.5.1;  and  some  constant  time  for  steps  4.4  and  4.5.2.  Thus 
algorithm  A4  requires  0(log(n+m))  time  for  detection. 
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'TaUe  I  St  wage  structure  of  clockwise  convex  chains  in  Fig.  7 
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*♦  INFO[i]  =  j,  i.e.  information  on  vertex  i  stored  at  location  j. 
*  1  means  a  branch  point  and  0  a  non-branch  point. 


As  for  algorithm  A5,  the  detection  time  complexity  is  0(n+m),  n+m 
being  the  total  number  of  vertices.  The  detection  step  is  then  of  the 
same  order  of  complexity  as  the  preprocessing  and  it  does  not  help. 

The  important  steps  of  edgorithms  A1  to  A4  are  summarized  in  the  fol¬ 
lowing  lemmas. 

Lemma  1 

The  sector  containing  a  vector  whose  origin  is  a  kernel  point  k  of  a 
preprocessed  star  polygon  can  be  found  in  O(logn). 

Lemma  2 

The  Intersection  of  a  line,  a  line  segment,  or  a  convex  chain  with  m 
vertices  and  a  convex  polygon  with  n  vertices  can  be  found  in  O(logn), 
O(logn)  and  0(log(n+m))  time  [ChaBO,  DobBl]. 

Lemma  3 

Clockwise  (counter-clockwise)  convex  chains  of  all  spans  of  a  star 
polygon  can  be  found  implicitly  in  0(n)  time  and  stored  implicitly  in 
0(n)  space. 


Lemma  4 
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The  visible  chain  of  a  convex  polygon  with  respect  to  a  point  k  owLtside 
ol  the  pc^ygon  can  be  found  in  0(logcnO,  where  m  is  the  number  oT  ver¬ 
tices  of  the  polygon  [Pre79]. 

LCTomah 

The  vertex  or  vertices  of  partition  of  a  visible  chain  of  a  convex 
polygon  with  respect  to  a  point  k  outside  of  the  polygon  and  a  refer¬ 
ence  line  can  be  found  in  O(logn)  time  using  Fibonacci  search  (see 
[ChaSO]  for  details),  where  n  is  the  number  of  vertices  of  the  polygon. 

5.1.  Critical  portion  of  a  convex  chain 

In  the  evaluation  of  the  time  complexity  above  for  algorithms  A3  and 
A4,  we  assume  that  we  can  do  binary  search  on  all  convex  chains.  This 
would  be  possible  only  if  we  are  willing  to  use  O(n^)  space  which,  in 
turn,  would  require  O(n^)  preprocessing  time. 

It  is  our  intention  now  to  prove  that  a  whole  convex  chain  would  not  be 
necessary.  In  fact,  the  useful  portion  of  a  convex  chain  from  a  given 
vertex  ends  at  the  first  branching  vertex  in  both  clockwise  and 
counter-clockwise  order  if  the  proper  choice  of  convex  chain  could  be 
found.  So  the  strategy  here  is  to  exploit  this  property  and  show  that 
such  a  chain  can  be  found  in  O(logn)  and  prove  that  it  is  sufficient. 
Before  doing  so,  we  state  the  problem  and  related  solutions  in  the  fol¬ 
lowing  lemmas. 

Lemmas 

The  detection  of  an  edge  intersection  of  a  star  polygon  S  and  a  line 
segment  Is  which  has  its  endpoints,  lep  and  rep,  outside  of  S  and  does 
not  intersect  the  convex  hull  of  S  is  completely  determined  from  the 
intersection  of  Is  and  2  oppositely  oriented  convex  chains  or  2  oppo¬ 
sitely  oriented  closest  convex  chains.  The  2  convex  chains  are:  the 
clockwise  convex  chain  ending  at  one  of  vertices  of  the  sec¬ 

tor  containing  lep  {vj-  if  line[rep,lep)  intersects  line[vj  and  Vj^-i 
otherwise)  and  starting  at  a  convex  hull  vertex  of  S;  and  the  counter- 
clookwise  convex  chain  starting  at  one  of  vertices  of  the  sec¬ 

tor  containing  rep  (1/^+1  if  line[lep,rep)  intersects  line[i'i+i,Ui)  and  Vi 
otherwise)  and  ending  at  the  other  convex  hull  vertex,  where 
(k,lep,rep)  is  a  clockwise  oriented  triangle  and  k  is  a  kernel  point  of  S. 
The  2  closest  convex  chains  are  determined  from  the  endpoints  of  Is  - 
the  closest  clockwise  convex  chain  from  lep  and  the  closest  counter¬ 
clockwise  convex  chain  from  rep.  The  closest  clockwise  convex  chain 
starts  at  the  nearest  branch  point  in  clockwise  order  from  lep  and  the 
closest  counter-clockwise  convex  chain  ends  at  the  nearest  branch 
point  in  counter-clockwise  order  from  rep. 

Lemma? 

The  detection  of  an  edge  intersection  of  a  star  and  a  visible  convex 
chain  contained  in  some  span  of  the  star  is  completely  determined  by 
the  intersection  of  the  visible  convex  chain  and  2  oppositely  oriented 
convex  chains  determined  by  the  sector  containing  the  vertex  of  par¬ 
tition  of  the  visible  convex  chain  or  2  oppositely  oriented  cl r>sp.9;t  con¬ 
vex  chains  determined  by  the  vertex  of  partition  of  the  visible  convex 
chain. 
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Tbe  case  in  which  there  are  2  vertices  of  partition  is  similarly  deter¬ 
mined  ELS  in  above  for  the  head  and  tail  portions  of  the  convex  chains 
and  from  Lemma  6  for  the  line  segment  portion. 

Lenmiafi 

The  planar  subdivision  induced  by  a  star  polygon  and  its  set  of  clock¬ 
wise  (counter-clockwise)  convex  chains  can  be  triangulated  uniquely 
in  0(n)  time. 

From  Obs.  3  we  know  that  the  interior  of  a  star  polygon  can  be  triangu¬ 
lated  in  0(n)  time. 

The  convex  cheiins  lay  the  ground  work  for  the  triangulation  of  the 
regions  inside  them  and  outside  of  the  boundciry  of  the  star  polygon. 
yfe  can  consider  triangulation  of  each  span  independently  of  the  other 
as  they  can  be  triangulated  in  a  similar  manner,  (Note  that  the  outer¬ 
most  convex  chain  is  an  edge.)  Beginning  with  the  vertex  which  is  the 
last  one  in  the  span  in  counter-clockwise  order  and  using  it  as  the 
anchored  point  of  tri angulation,  connect  an  edge  to  each  vertex  of  the 
clockwise  convex  chain  visible  on  its  right  side  that  is  currently  not 
connected  to  it.  (The  subdivision  induced  by  counter-clockwise  con¬ 
vex  chains  is  similarly  triangulated  beginning  with  the  first  vertex  in 
the  span  in  the  counter-clockwise  order.)  Since  each  vertex  is  visited 
at  most  twice  and  the  convex  chains  Eire  already  avedlable,  it  then 
takes  0(n)  time  to  do  the  required  triangulation.  (See  Rig.  8  for  the 
result  of  the  trieingulation  of  the  star  polygon  and  its  clockwise  convex 
chains  in  Fig.  7.) 

Lemma  9 

The  triangular  subdivision  problem  can  be  solved,  with  0(n)  space  Euid 
0(n)  preprocessing  time,  in  O(logn)  time  [KirSl]. 

From  Lemma  3,  we  know  how  to  find  the  convex  chains  implicitly  in 
0(n)  time  and  use  0(n)  space  by  storing  only  each  of  the  vertices  once 
and  only  branching  vertices  for  overlapping  convex  chains.  The 
remaining  problem  of  now  to  find  the  closest  convex  chain  from  a 
given  point  in  O(logn)  time  can  be  solved  using  Lemmas  8  and  9. 
Lemma  8  guarEintees  that  only  an  extra  0(n)  preprocessing  time  is 
necessary  to  solve  triangular  subdivision  problem  of  planar  subdivi¬ 
sion  induced  by  a  star  polygon  and  its  convex  chains.  Lemma  9 
guarantees  us  that  we  can  find  the  closest  convex  chain  in  O(logn) 
time,  thus  obtaining  the  desired  detection  time  complexity  without 
having  to  pay  the  price  of  O(n^)  preprocessing  time  and  space. 

6.  Conclusion 

The  results  of  this  paper  are  summarized  in  Table  11  together  with 
some  other  known  relevEint  results.  In  the  algorithms  presented 
above,  specifically  in  the  preprocessing  step,  we  use  some  of  the 
known  results  on  how  to  find  kernel  and  convex  hull  and  how  to  obtain 
optimal  searching  in  a  triangular  subdivision.  In  using  them,  one 
must  be  careful  that  the  Input/output  formats  of  the  known  algo¬ 
rithms  are  compatible  or  else  some  additional  order  of  time  may  be 
needed  to  convert  from  one  standard  representation  to  the  other.  It  is 
therefore  necessary  to  start  with  the  proper  representation  to  avoid 
extra  conversion  cost.  In  our  case,  the  representation  must  be  compa- 
tiMe  with  the  edge-ordered  represention  in  [Kirei]. 


-24- 


Fig.  6  Result  of  triangulating  the  star  polygon  and  its  clockwise  con¬ 
vex  chains  in  Fig.  7. 


Note  that  though  algorithm  A3  could  be  used  to  find  intersection 
between  a  convex  polygon  and  a  star  polygon,  careful  analysis  indi¬ 
cates  that  O(nlogn)  time  might  be  necessary  in  the  worst  case,  where 
both  polygons  have  0(n)  vertices. 

In  this  paper,  we  show  how  to  exploit  the  visibility  property  of  a  star 
polygon  and  its  convex  chains  to  facilitate  intersection  detection.  We 
must  note  here  that  two  edges  or  two  sequences  of  edges  could  be 
hanging  onto  a  convex  chain  where  intersection  test  is  being  made. 
However,  only  when  a  convex  chain  edge  intersection  is  detected 
should  the  proper  sequence  of  star  polygon  edges  be  tested.  (See  Fig. 
9.)  This  idea  could  be  exploited  for  other  structurally  well-defined 
polygons  and  maybe  also  for  3-dimensional  objects. 
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Table  n  Intersection  detection  complexity 


obiect 

star  polygon 

cony  ex  polygon 

point 

Odogn) 

Sha70l 

O(logn) 

r  Sha7Bl 

line 

O(logn') 

O(logn) 

[ChaBOl 

line  segment 

O(logn) 

Oflogn)  [ChaBOl 

convex  polygon 

O(logn) 

M 

Oflogn)  [ChaBOl 

star  polygon 

0(n^ 

O(logn) 

L^] _ 

♦  These  results  are  shown  in  this  paper. 


fig  9  Some  proper  sequences  of  a  star  polygon’s  edges  associated 
with  vertices  of  convex  chain(t'i,i>5,x'iQ,i^  12). 
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